LO Phonon-Induced Exciton Dephasing in Quantum Dots: An Exactly Solvable Model 
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It is widely believed that, due to its discrete nature, excitonic states in a quantum dot coupled 
to dispersionless LO phonons form everlasting mixed states (exciton polarons) showing no line 
broadening in the spectrum. This is indeed true if the model is restricted to a limited number of 
excitonic states in a quantum dot. We show, however, that extending the model to a large number 
of states results in LO phonon-induced spectral broadening and complete decoherence of the optical 
response. 
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Among different mechanisms of the optical decoher- 
ence in semiconductor quantum dots (QDs), interaction 
between excitons and lattice vibrations (phonons) is the 
most important one at low carrier densities, leading to 
a temperature-dependent linewidth in optical spectra. 
The discrete nature of excitonic levels in a QD gave rise 
to the so-called phonon bottleneck problem 0]: When 
the phonon energy does not fit the level separation, real 
phonon-assisted transitions between different excitonic 
states are not allowed. On the other hand, the mea- 
sured optical polarization in QDs shows a partial ini- 
tial decoherence and a temperature-dependent exponen- 
tial decay at larger times |2|. The bottleneck problem 
is partly removed for acoustic phonons as they have an 
energy dispersion and thus contribute to carrier relax- 
ation Moreover, even apart from the resonance be- 
tween the phonon energy and the level separation, acous- 
tic phonons, due to their dispersion, are shown to be 
responsible for pure dephasing induced by virtual pro- 
cesses |4|. 

Longitudinal optical (LO) phonons, in turn, have a 
negligible dispersion which leads to a quite different be- 
havior: Excitons in QDs and LO phonons are always in a 
strong coupling regime in the sense that they form ev- 
erlasting polarons with no spectral broadening 
This also makes any approximate perturbative approach 
to the exciton-LO phonon problem inappropriate. In- 
deed, for a few excitonic levels in a QD coupled to bulk 
LO phonons, a self-consistent second Born approxima- 
tion for the self energy developed in [ToL Hl| results in 
a line broadening which is fully artificial. This artefact 
is refuted by the exact solution of this problem 0, 
showing that the spectrum consists exclusively of dis- 
crete unbroadened lines. Also, a Gauss spectral lineshape 
was found in the quadratic coupling model by truncating 
the cumulant expansion in second order |12( . Again, the 
exact solution of this model shows no spectral broaden- 
ing [13. 

In this Letter we show that pronounced LO phonon- 
induced dephasing and spectral broadening do neverthe- 
less exist in QDs. This broadening is calculated micro- 



scopically by inclusion of infinitely many excitonic states 
in a QD that has never been done before. Thus a quali- 
tatively new source of the dephasing in QDs is found. 

The full problem of excitons in a QD linearly coupled 
to the LO phonon displacement can be solved exactly 
only for a very limited number of excitonic states |9(. 
The major obstacle are off-diagonal terms in the exciton- 
phonon interaction which couple different excitonic states 
in a QD. In order to take into account as many excitonic 
states as we like, we have derived microscopically an ef- 
fective exciton-phonon Hamiltonian which has only level- 
diagonal terms and thus allows an exact solution for an 
arbitrary number of states. This effective Hamiltonian, 
however, preserves the main features of the original prob- 
lem, since the off-diagonal terms are also intrinsically rep- 
resented: By means of a unitary transformation they are 
mapped into diagonal terms giving rise to an exciton- 
phonon cou pling quadratic in the phonon displacement 
operators jj,|lj|. 

Concentrating on the ground exciton state |1), the ef- 
fective Hamiltonian takes the form (H = 1): 

H = u; J2a{a<i+(Ei + V L + V Q )\l}(l\, (1) 
q 

Vl = X>"(q)K + «-q)> (2) 
q 

V Q = -^Q(p,q)(a P + a t _p)K + « t _ q ),(3) 
pq 

Q(p,q) = 2^ En ~^ , M lw (p)Af nl (q) 

= ^F n (p)F„(q), (4) 

where u>o is the dispersionless LO phonon frequency 
and E n is the bare transition energy of a single-exciton 
state \n) in a QD. The kernel Q of the quadratic cou- 
pling to LO phonons has the same form as the ef- 
fective scattering matrix used in Ref. Q to describe 
the phonon modes bound to a neutral donor. It 
is derived from the levcl-nondiagonal matrix elements 
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Mi„(q) oc q _1 (l|e 4ql " e — e lqr ''|n) of the linear exciton- 
phonon (Frohlich) interaction, up to second order. Us- 
ing its factorization property, Q is expressed in Eq. Q 
in terms of functions F n (q). A quadratic coupling 
model similar to Eqs. I^-^J has been already successfully 
used for calculation of the exciton dephasing in InGaAs 
QDs, induced by acoustic phonon-assisted virtual transi- 
tions 0| . 

The method developed in allows us to find the exact 
solution of the Hamiltonian using the cumulant 

expansion. The linear optical polarization has the form 

P(t) = 0(t) exp[-iE 1 1 + K L (t) + K Q {t) + K M {t)\ , (5) 

where the linear and quadratic cumulants are 

K L (t) = [(2N+l)(cosu} t-l)+i(ujQt-smu; t)} 

x^|Af n (q)/ W o| 2 , (6) 
q 

K Q® = -^EHl-iA^W]. (7) 

Here A„ are the eigenvalues of the matrix 

A nm = ^2F n (-q)F m (q), (8) 
q 

while Xj (t) are the eigenvalues of the Fredholm problem 

[ dT'D{T-T l )u j {T l ;t) = X j {t)u J (T\t). (9) 
Jo 

Equation © can be solved numerically, as it has been 
done for acoustic phonons 0. However, for dispersion- 
less optical phonons the propagator D is q-independent: 
D(t) = (N+l)e-^°l T l + 7V e ^oM [ w here = l/(e f3u °-l) 
and j3 = 1/ksT] and thus allows an analytic solution 
of Eq. |J3J which can be found in Ref. ^3|. Finally, the 
mixed cumulant ATm(^) in Eq. (JSJ is found in the same 
way as the quadratic one but additionally requires the 
eigenvectors of the matrix A nm and the eigenfunctions 
Uj of Eq. ©. It turns out, however, that it leads to only 
tiny corrections: \Km/Kq\ < 10~ 7 in CdSe QDs consid- 
ered here, and thus can be safely neglected. 

We have already used the model Eqs. and its 

solution Eqs. El-EH in case of two excitonic levels in an 
InGaAs QD [ljj to show that a truncation of an infinite 
series in the cumulant as done in Ref. ^3] leads to an ar- 
tificial Gauss decay of the optical polarization. Only tak- 
ing into account an infinite number of all-order diagrams 
gives the correct result which is an almost perfectly peri- 
odic time-dependent polarization. In the present Letter, 
using the same model, we include now (infinitely) many 
excitonic states in a QD and show that this leads to a 
qualitatively new effect: polarization decay and spectral 
broadening. 



This effect takes place in any type of semiconductor 
QDs. For the sake of illustration, we concentrate on the 
polar material CdSe and a simple model of a QD. Mate- 
rial parameters (Frohlich coupling constants a e = 0.47, 
ah = 0.88; Hujq = 24meV) and the QD model are the 
same as in Ref. @, but the parabolic confinement po- 
tentials are taken as isotropic. For the Gauss localiza- 
tion length of electron and hole we use here three sets of 
parameters: (A) l e = lh = 2.5 nm; (B) l e = 2.71 nm, 
l h = 2.2 nm (as in Ref. 0); and (C) l e = 4.14 nm, 
lh = 3.36 nm, which gives half the level distance com- 
pared to case (B). We decided to neglect the Coulomb 
interaction as it would only renormalize the transition 
energies E n and induce minor changes in the polariza- 
tion dynamics of the ground exciton state. 

The amplitude of the time-dependent linear polariza- 
tion after impulsive excitation in a CdSe QD at T = 
150 K is shown in Fig. 1 for set (A). In this case, due 
to the full charge neutrality of the electron-hole pair, 
Afn(q) = holds, and the linear term Vl vanishes. 
Thus we concentrate at the moment on the quadratic 
term Vq , Eq. © , which is in fact the only source of de- 
phasing. For two excitonic levels [i.e., only one excited 
state n = 2 contributes to Eq. J3J|, and the matrix A nm , 
Eq. |JSJ|, reduces to a scalar] the polarization has only tiny 
oscillations close to unity. For five levels, it decays and 
then revives almost to unity after 20 ps. For 30 levels, 
the amplitude decays already dramatically and then os- 
cillates irregularly around a small value. Such a behav- 
ior has been already discussed before and called quasi- 
dephasing 0, 0] . Finally, if we take into account all ex- 
citonic states, the polarization decays strictly: At 100 ps 
its amplitude drops to 10 -30 (but is shown in Fig. 1 only 
up to 10~ 8 ). In practice, of course, we truncate the infi- 
nite matrix A nm and check that a larger matrix produces 
no changes within the time interval considered. 

To understand better why the optical polarization be- 
comes decaying when more and more exciton (electron- 
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FIG. 1: Amplitude of the linear optical polarization of the 
ground exciton state at T — 150 K, taking into account two, 
five, thirty, and all exciton levels in a CdSe QD (set A). 
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hole pair) states are taken into account, we calculate 
the linear absorption spectrum, i.e. the real part of the 
Fourier transform of P{t). It is hopeless, however, to find 
the Fourier transform numerically if the function does not 
decay at all or decays very weakly. That is why we use 
here a different method: We diagonalize exactly the ex- 
cited state Hamiltonian Hi which is related to the full 
one, Eq. fl), as H = #o|0)(0| + #i|l)(l|. Such a diago- 
nalization is very easy if Vl = 0. In case of two levels Hi 
becomes 



H^ lev = nB^B + Ei+cj ^2b%, (10) 

n = ^ 2 -2 Wo ^F 2 (-q)^ 2 (q), (11) 
B = ^FaCcO^+Oq + ^al,) (12) 



where 



are, respectively, the new phonon frequency and annihi- 
lation operator, £± = (fl±uj )/2^/uJo^J2 cl l-^Mq)! 2 , and 



Ei is the polaron shifted transition energy. 

Note that out of a continuum of degenerate LO modes 
only a single phonon mode couples to the QD and pro- 
duces a new, bound mode |l ll B with frequency Q, while 
all the other modes b^, which show up in the last term of 
Eq. l|lUjl. are decoupled. They are simply some orthogo- 
nal linear combinations of the former phonon modes a q 
and have the same old frequency u>q. The linear polar- 



ization P(t) = 
way as in Ref. 



0(t)(e iHot e~ iHlt ), calculated in the same 



then takes the form: 



'(*) = m 



otnm exp{i(nwo — mfi)i}, (13) 



where a r 



(1 



e _/3n " |(n|m)| 2 , and (n|m) are 
the projections of new phonon states |m) into old ones 
n) (to be distinguished from exciton states \n)), which 
are calculated recursively. 

Generalization to TV excitonic levels is straightforward: 
Hi is now diagonalized to TV — 1 new phonon modes: 



2cj n A 



(14) 



with A„ being the eigenvalues of the (TV — l)-dimensional 
matrix Eq. JSJ). In accordance with Eqs. (JSJ and @, the 
full linear polarization can be written as a product 



P(t) = ]JPAt) 



(15) 



of functions Pu{t) due to each individual phonon mode 
given by the same Eq. I|13fl , where f2 and a nm are re- 
placed, respectively, by Q, u and o^ m . 




-1.5 -1.0 

co/ox 



-1.5 -1.0 

(0/(0„ 



FIG. 2: Evolution of the fine structure in the absorption spec- 
trum of a QD by accounting for more and more excitonic levels 
[linear scale, set (A), E\ is taken as zero of energy]. 



For two levels, the absorption spectrum, already dis- 
cussed in Ref. [l3|, represents a set of discrete lines. 
There are two-phonon satellites around the zero-phonon 
line (the standard one-phonon satellites are absent here 
as Vl — 0), but more important is the fine structure 
shown in Fig. 2 (a). This fine structure is due to the dif- 
ference between the old and the new phonon frequency 
loq — Q and manifests itself in the time domain as rapid 
oscillations with the period of 0.5 ps (Fig. 1). Including 
the third excitonic level brings in an additional frequency 
and more lines in the spectrum, Fig. 2 (b). With five lev- 
els (and four new phonon frequencies), there is already a 
plenty of closely lying delta lines, Fig. 2 (c). Finally, if we 
include all exciton levels, these lines merge and produce 
a continuous broadening, Fig. 2 (d). In the time domain, 
the polarization represents a superposition of individual 
oscillations [see Eqs. I|15f) and (|13fl ] which interfere de- 
structively due to incommensurability of their frequen- 
cies and thus lead to dephasing. 

In contrast to the present model of a QD with discrete 
exciton levels only, in realistic QD systems there is also a 
wetting-layer continuum. This continuum could seem to 
be a more probable candidate for producing the spectral 
broadening [l6|. It turns out, however, that in typical 
QDs, there are usually enough discrete exciton levels to 
provide an appreciable polarization decay. In fact, taking 
into account 30 levels is already sufficient to have the 
dephasing time very close to the exact one (i.e., with 
all levels, see Fig. 1), but at the same time it requires 
the heterostructure potential band offset to be less than 
1.0 eV for the electron and 0.5 eV for the hole 01- We 
account for all discrete levels in the harmonic potentials 
just for consistency, in order to show that a complete 
decay of the polarization is indeed achieved. 
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FIG. 3: Absorption spectrum of a CdSe QD at T = 300 K, cal- 
culated for two different types of the bare-exciton level struc- 
ture [sets (B) and (C), see text]. Bound phonon frequencies 
Qv and levels E n of exciton excited (ground) states are shown 
in the insets by thin (thick) horizontal lines. 



The full absorption spectrum of a CdSe QD is shown 
in Fig. 3. As we have in sets (B) and (C) different Gauss 
lengths, l e ^ lh, the linear coupling Vl contributes too, 
which gives surprisingly only a slight modification in the 
spectrum. For set (B) all the excited levels are well above 
Hloq. All the bound-phonon frequencies Sl v are condensed 
just below ojq, except one which is well separated and is 
in fact responsible for the spectral fine structure. This 
fine structure is absent if we take a larger dot with closer 
levels (Fig. 3, red curve and inset C). Here, one excited 
level is found below Hloq, and even two bound phonons 
are well resolved and wash out the fine structure. Thus, 
changes in the exciton level structure modify the spec- 
trum but obviously do not affect our principal conclusion 
on the finite spectral width, and even do not alter much 
the linewidth itself. 
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FIG. 4: Spectral linewidth of the ground exciton state in a 
CdSe QD, set (B), as a function of temperature (solid curve) 
and the exponential fit with C = 20meV (dashed curve). 
Inset: Polarization amplitude at different temperatures. 



Since the decay is not strictly exponential (see the inset 
in Fig. 4), we have extracted the linewidth directly from 
the absorption and plotted in Fig. 4 the full width at half 
maximum (FWHM) as a function of temperature. At 
high temperatures, the FWHM varies from few to ten 
meV and can be well approximated by an exponential 
law with hajQ being the activation energy. Below 50 K 
the linewidth is practically unaffected by LO phonons, 
and the dephasing in QDs is determined by some other 
mechanisms like coupling to acoustic phonons , phonon 
anharmonicity |18|, and radiative decay. 

In conclusion, we have calculated the time-dependent 
optical polarization and the absorption spectrum of a 
quantum dot using the exactly solvable model of excitons 
quadratically coupled to LO phonons. We have extended 
this model to an arbitrary number of excitonic states in a 
quantum dot and (using its exact solution) demonstrated 
that taking into account infinitely many states results in 
an LO phonon-induced exciton dephasing and spectral 
broadening. 
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